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PREFACE 

This  document  presents  a  general  mathematical  development  of  the 
theory  of  vapor  contamination  hazards  to  nuclear,  biological  and  chemical 
(NBC)  collective  protection  systems.  An  important  general  result  is 
derived  and  then  illustrated  by  applying  it  to  selected  examples  of 
special  interest  in  studies  of  NBC  defense. 

This  work  was  performed  under  the  Naval  Aircraft  Combat  Survivability 
Program's  Task  F-82-11  "Aircraft  Protection/Decontamination  Requirements 
for  Naval  Air  Operations  in  a  CBR  Environment." 
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GLOSSARY  OF  SYMBOLS  AND  ACRONYMS 


A  fixed  control  volume  or  region ,  e.g.,  the  Interior  of  a 
collective  protection  system 

3 

Volume  in  m  of  R 

3 

Instantaneous  volume  rate  of  flow  in  m  /min  of  inflow 
stream  i 

3 

Instantaneous  volume  rate  of  flow  in  m  /min  of  the  outflow 
stream;  see  Eq.  (2) 

Number  of  inflow  streams 

3 

Instantaneous  concentration  in  mg/m  of  contaminant  vapor 
in  R;  see  Eq.  (5) ,  (9)  and  (11) 

Instantaneous  mass  in  mg  of  contaminant  vapor  in  R; 
see  Eq.  (1)  and  (4) 

Instantaneous  stopping  power  or  efficiency  of  filter/ 
absorber  unit  j;  see  Eq.  (3) 

3 

Instantaneous  concentration  in  mg/m  of  contaminant 
vapor  at  the  output  of  filter /absorber  unit  j 

3 

Instantaneous  concentration  in  mg/m  of  contaminant 
vapor  in  inflow  stream  i 

3 

Instantaneous  volume  rate  of  flow  in  m  /min  through 
filter/absorber  unit  j 

Number  of  filter/absorber  units 

Number  of  evaporating  sources 

Instantaneous  rate  of  evaporation  in  mg/min  of 
contaminant  vapor  source  k  at  time  t 

Weighted  average  inflow  stream  contaminant  vapor 
concentration  in  mg/m  (see  Eq.  (6.1) 

3 

Total  rate  of  flow  in  m  /min  through  the  filter/absorber 
units;  see  Eq.  (6.2) 
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Weighted  average  filter /absorber  stopping  power  or 
efficiency;  see  Eq.  (6.3) 

Total  rate  of  evaporation  in  mg/min  of  the  contaminant 
vapor  sources  in  R;  see  Eq.  (6.4) 

Instantaneous  fractional  ventilation  rate;  see  Eq.  (8.1) 

Instantaneous  rate  of  introduction  of  contaminant  vapor 
per  unit  volume  of  R  in  mg  *  min  •  m  ;  see  Eq.  (8.2) 

A  time  to  which  initial  concentrations  and  dosages  are 
referenced 

3 

Vapor  dosage  in  mg-min/m  accumulated  in  R  from  time  t 
to  time  t;  see  Eq.  (10)  ° 

3 

Total  vapor  dosage  in  mg-min/m  accumulated  in  R  over 
all  time;  the  limit  of  Dp(t)  as  t  -*•  °°;  see  Eq.  (13) 
and  (15) 

The  limit  of  CR(t)  as  t  +  “ 

A  time  chosen  large  enough  that  B(z)  <  ea  for  z  _>  T 
The  time  at  which  a  change  in  the  value  of  a  occurs 
Value  of  a  for  t  <  t  <  t. 

Value  of  a  for  t,  *-  t  <  " 
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INTRODUCTION 

Our  purpose  is  to  demonstrate  a  useful  and  interesting  general 
result  in  the  theory  of  vapor*  contamination  hazards  to  collective 
protection  systems.  The  combination  of  the  generality  of  the  result  and 
the  simplicity  of  its  proof  may  be  new — at  least  I  do  not  recall  having 
seen  either  of  them  before.  In  any  event,  both  the  general  result  and 
its  elegant  proof  deserve  to  be  more  widely  known.  It  can  easily  be 
used  to  obtain  for  particular  cases  the  following  striking  results. 

Example  1. — Suppose  that  shortly  before  an  NBC  attack  a  collective 
protection  shelter  with  no  filtered  air  supply  is  imperfectly  sealed  so 
that  a  small  but  persistent  flow  of  outside  air  enters  the  shelter.  If 
this  sealed  posture  of  the  shelter  is  maintained  indefinitely,  then  the 
total  vapor  dosage  inside  the  shelter  is  the  same  as  the  total  vapor 
dosage  outside  the  shelter.  Furthermore,  this  result  holds  irrespec¬ 
tive  of : 


1.  The  size  of  the  shelter  (provided  its  volume  is  constant) 

2.  The  leakage  flow  rate  (provided  it  is  a  positive  constant) 

3.  Changes  with  time  of  the  contamination  vapor  density  in 
the  outside  air  (provided  it  approaches  zero  as  time  passes) . 

Example  2. — Suppose  that  some  liquid  contaminant  is  inadvertently 
brought  into  a  previously  uncontaminated  collective  protection  shelter 
that  is  provided  with  a  perfectly  filtered  air  intake  and  has  no  leaks 
from  the  outside  air.  Then  the  total  vapor  dosage  in  mg-min  per  cubic 
meter  inside  the  shelter  is  equal  to  the  total  mass  in  mg  of  contaminant 
that  eventually  vaporizes  divided  by  the  shelter's  volume  ventilation 
rate  in  cubic  meters  per  minute.  Furthermore,  this  result  holds  irrespec¬ 
tive  of: 


*For  brevity,  all  airborne  contaminant  (whether  vapor  or  aerosol)  is 
referred  to  as  "vapor."  Similarly,  "vaporizing"  includes  re-aerosolizing 
and  so  forth. 
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1.  The  size  of  the  shelter  (provided  its  volume  is  constant) 

2.  The  volume  ventilation  flow  rate  (provided  it  is  a  positive 

constant) 


3.  Changes  with  time  of  the  liquid  contaminant's  vaporization 
rate  (provided  it  approaches  zero  as  time  passes). 

The  generalization  of  these  particular  cases,  and  the  result  which 
will  be  demonstrated  in  the  following  sections,  can  be  stated  in  words 
as  follows.  Consider  a  collective  protection  shelter  whose  volume  is 
fixed  and  whose  volume  ventilation  rate  in  cubic  meters  per  minute  is 
kept  constant.  Then  the  total  vapor  dosage  in  mg-min  per  cubic  meter 
inside  the  shelter  is  equal  to  the  total  mass  of  contaminant  vapor  in  mg 
that  is  initially  within  or  subsequently  introduced  into  the  shelter, 
divided  by  the  shelter's  volume  ventilation  rate.  Furthermore,  this 
result  holds  irrespective  of: 

1.  The  size  of  the  shelter  (provided  its  volume  is  constant) 

2.  The  shelter's  volume  ventilation  flow  rate  in  cubic  meters 
per  minute  (provided  it  is  a  positive  constant) 

3.  Changes  with  time  of  the  rate  in  mg  per  minute  at  which 
contaminant  vapor  is  introduced  into  the  shelter  by  vaporization  of 
internal  sources,  by  the  influx  of  imperfectly  filtered  outside  air,  or 
both  (provided  that  the  total  rate  of  introduction  of  contaminant  vapor 
approaches  zero  as  time  passes) 

4.  The  initial  concentration  of  contaminant  vapor  within  the 

shelter. 

The  presentation  of  this  general  result  begins  with  a  description 
of  the  even  more  general  case  in  which  the  shelter's  ventilation  rate 
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DESCRIPTION  OF  THE  GENERAL  SITUATION 


Consider  the  region  R  shown  in  Figure  1.  R  could  be  any  fixed 

3 

control  volume  or  region  with  a  constant  volume  of  V  m  ,  but  it  usually 

is  interpreted  to  be  the  interior  of  a  collective  protection  system. 

Assume  that  any  contaminant  vapor  in  R  is  well-mixed  so  that  its  vapor 

3 

concentration  CR(t)  in  mg/m  is  the  same  at  all  points  of  R.  Thus,  if 
QR(t)  is  the  mass  in  mg  of  the  contaminant  vapor  in  R  then 

QR(t)  =  V  CR(t)*  (1) 

If  the  contaminant  vapor  is  not  well-mixed  there  may  be  "hot  spots"  of 
higher  than  average  concentration  in  R  which  are  balanced  by  "cold 
spots"  of  lower  than  average  concentration. 


Suppose  that  R  contains  K  pieces  of  solid  or  liquid  contaminant  and 
that  piece  k  vaporizes  at  the  rate  GR(t)  mg/min.  Suppose  that  R  is  fed 

by  I  influx  streams  and  that  influx  stream  i  has  a  volume  flow  rate  of 

3  3 

F^t)  m  /min  and  a  contaminant  vapor  concentration  of  C^(t)  mg/m  . 

Then  the  total  rate  of  volume  flow  into  R  through  the  influx  streams  is 

I 

F  (t)  =  I  F  (t),  (2) 

i=l  1 


Suppose  also  that  R  is  equipped  with  J  filter/absorber  units  and  that 

3 

unit  j  has  a  volume  flow  rate  of  F  . (t)  m  /min  and  that  its  output 

Aj  ^ 

stream  has  a  contaminant  vapor  concentration  of  C^(t)  mg/m  .  Assume 
that  filter /absorber  unit  j  has  a  stopping  power  or  efficiency  ratio 


E^t)  -  [CR(t)  -  CAj(t)]/CR(t)  (3) 

that  is  a  function  of  time  only.  Assume  that  there  is  no  appreciable 
gain  or  loss  of  volume  by  filtration  or  by  contaminant  vaporization,  and 
let  all  flows  be  practically  incompressible. 
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Then  the  ness  of  contaminant  vapor  in  R  at  time  t  must  satisfy  the 
mass  balance  equation 


K  I 

Qp  (t+dt)  -  Qp(t)  +  I  G.  (t)  dt  +  E  C.Ct)  F,(t)  dt 
R  R  k-1  k  i-1  1  1 


J 

-  Z 

j-1 


V‘>  [CRCt)  -  CAJ(e)1  dt 


-  F0(t)  CR(t)  dt.  (4) 

By  using  equations  (1) ,  (2)  and  (3) ,  equation  (4)  may  be  written  as 
d  C_(t) 

V  -JJ - Fo(t)  C0(t)  -  FAo(t)  CR(t)  E(t) 

-  F0(t)  CR(t)  +  G0(t)  (5) 


where  we  have  put 


co(t) 

FAo(t> 

E(t)  - 

Go(t) 


Z  F.Ct)  C . (t) /F  (t) 
i-1  1  1  ° 


2  FA,(t) 


j-1 


Aj 


Z  FAA(t)  E, (t)/F._(t) 


j-1 


Aj  j  Ao 


E  Gk(t>’ 

k-1  K 


(6.1) 


(6.2) 


(6.3) 

(6.4) 


Note  that  E(t)  is  dimensionless  and  represents  a  sort  of  weighted  average 
filter/absorber  efficiency  factor,  and  C  (t)  is  a  weighted  average  inflow 

3  ° 

vapor  concentration  expressed  in  mg/m  .  Equation  (5)  can  be  written  as 
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d  C^t) 

— "  -a(t)  CR(t)  +  8(t)  (7) 

where 

a(t)  -  [Fq (t)  +  FA<j(t)  E(t)]/V  (8.1) 

and 

6(t)  -  [GQ(t)  +  FQ(t)  CQ(t)J/V  (8.2) 

where  a(t)  and  g(t)  are  functions  of  t  but  do  not  depend  explicitly  on 
other  quantities.  Physically,  ct(t)  and  g(t)  are  non-negative.  We  call 
«(t)  the  instantaneous  fractional  ventilation  rate,  and  V  a(t)  the 
instantaneous  volume  ventilation  rate  of  the  region  R  (or  of  the  collec¬ 
tive  protection  system  which  it  represents).  Observe  that  g(t)  is  the 

instantaneous  rate  of  introduction  of  contaminant  vapor  per  unit  volume 

-1  -3 

of  R  and  has  the  dimensions  mg  •  min  •  m  .  V  B(t)  is  the  Instantaneous 
rate  in  mg/min  at  which  contaminant  vapor  is  introduced  into  R. 


Equation  (7)  is  a  linear  differential  equation.  Its  general  solution 
is  presented  in  many  standard  texts  on  ordinary  differential  equations 
and  may  be  written  as 


-  /t  a(u)du 
o 


CR(t) 


E 


f t  a(s)ds 
o 


Cp (t  )  +  f  B(z)  e 
K  o  t 

o 


dz  !  (9) 


where  C_(t)  is  the  concentration  in  R  at  time  t  . 

K  O  3 

dosage  or  cumulative  C-t  product  in  R  in  mg-min/in 


By  definition,  the 
after  time  t  is 


t 

DR(t)  -  /  CR(y)dy.  (10) 

Co 

At  this  level  of  generality  we  have  not  yet  been  able  to  deduce  any 
particularly  noteworthy  properties  of  the  total  vapor  dosage.  Accordingly, 
we  consider  instead  the  constant  ventilation  rate  situation. 
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GENERAL  RESULT  FOR  CONSTANT  VENTILATION  RATES 

When  the  ventilation  rate  a  is  constant  Eq.  (9)  reduces  to 

t 

CR(t)  -  CR(tQ)exp(-  a(t-tQ))  +  f  8(z)exp(-  a(t-z))dz  (11) 

t 

o 

Our  main  object  is  to  demonstrate  the  following  general  result  for 
collective  protection  shelters  with  constant  positive  ventilation  rates. 

THEOREM:  Suppose  a  is  a  positive  constant  and  '  (z)  is  such  that 


C_(»)  -  lim  C  (t)  «  0. 

t  -*•  «R 


(12) 


Then 


Dr(~)  -  lim  DR(t)  -  a 
t  -*■  0» 


-1 


[CR(t0)  +  f  6 (y)dyj . 


(13) 


PROOF:  When  a  is  a  positive  constant  we  may  solve  Eq.  (7)  for 
C„(y)  and  substitute  the  result  into  Eq.  (10)  to  obtain 


-1 


DR(t)  -<*  f t  [6(y)dy  -  dCR(y)J 


a_1[CR(to)  +  /  S (y)dy  -  C  (t)]. 
Lo 


(14) 


But  CR(“)  *  0  by  assumption,  and  Eq.  (13)  follows  immediately.  QED. 


Substituting  from  Eqs.  (8.1)  and  (8.2)  into  Eq.  (13)  we  may  write 


Dr(»)  -  [VCR(tQ)  +  £  (G(y)  +  FQ(y)  CQ(y))dy]  /  tFQ  +  FAE], 


(15) 
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where  by  the  assumptions  of  the  Theorem  the  denominator  is  a  positive 
constant  (i.e.,  independent  of  t).  Observe  that  by  Eq.  (1)  the  first 
term  of  the  numerator  is  the  initial  mass  of  contaminant  vapor  in  R. 

The  second  term  of  the  numerator  is  the  mass  of  contaminant  vapor 
subsequently  introduced  into  R  by  all  causes.  The  denominator  is  the 
volume  ventilation  rate  of  R.  Thus  the  conclusion  of  the  Theorem  can  be 
stated  in  words  as  given  in  the- Introduction. 

It  can  be  inconvenient  to  verify  that  B(z)  is  such  as  to  satisfy 
condition  (12).  The  following  Lemma  is  often  easier  to  apply. 

LHiMA:  Suppose  ct  is  a  positive  constant  and 

6(»)  =  lim  g(t)  *  0. 
t  00 

Then 

C_(»)  =  lim  Cw(t.)  -  0. 


PROOF:  Given  any  e  >  0  fix  T  so  large  that 


8(z)  <  za  for  z  _>  T. 

(Since  by  assumption  S(z)  ■+■  0  as  t  +  ®,  an  appropriate  T  can  always  be 

found).  Then  by.  breaking  the  integral  in  Eq.  (11)  at  t-T  we  see  that 


T  oi(z-to)  _  a(t_t  ) 

CR(t)  <  z [ 1-e”  a(t“T)]  +  [CR(to)  +  /  8 (z)e  dz]e  °  . 


Letting  t  -*■  00  in  this  last  formula  yields 


CR(°°)  <  c. 


But  since  z  may  be  any  positive  number,  we  must  in  fact  have 


CR(»)  =■  0. 


QED 


17 


NADC-82255-20 


DISCUSSION 

It  is  true  that  6(<»)  -  0  is  not  a  sufficient  condition  for  the 

convergence  of  t 

/  8(y)<ly 
t 

o 

to  a  finite  limit  as  t  +  «.  However,  since  8(y)  is  non-negative,  this 
integral  will  either  converge  to  a  finite  limit,  or  else  diverge  properly 
to  +  oo.  The  latter  situation  corresponds  to  the  eventual  introduction 
as  time  passes  of  an  infinite  mass  of  contaminant  vapor  into  R.  Eq.  (13) 
should  be  interpreted  as  being  a  relation  between  finite  quantities  when 
the  integral  on  the  right  converges,  and  as  asserting  that  when  the 
integral  diverges  properly  to  +  «>  so  does  DR(t). 

The  assumption  that  a  is  positive  is  essential  in  both  the  Lemma 
and  Theorem.  When  «  ■  0  Eq.  (11)  shows  that 

t 

CR(t)  -  CR(tQ)  +  f  B(z)dz, 

Co 

and  so  C  ( “)  need  not  vanish  even  when  6(z)  is  identically  equal  to 
K 

zero,  which  is  incompatible  with  the  conclusion  of  the  Lemma.  Sub¬ 
stituting  this  expression  for  CR(t)  into  Eq.  (10)  yields 

t  y 

DR(t)  -  CR(tQ)(t-to)  +  /  dy  /  6(z)dz, 

Co  Co 

which  shows  that  when  a  *  0  D„(t)  may  diverge  to  00  as  t  **■  00  even  when 

K 

8  (z)  is  identically  equal  to  zero,  which  is  incompatible  with  the 
conclusion  of  the  Theorem  as  expressed  by  Eq.  (13). 

The  assumption  that  ct  is  constant  is  also  essential  to  the  validity 

of  the  Theorem.  For  suppose  that  B(z)  is  identically  equal  to  zero, 

that  C-,(t  )  >  0,  and  that  a(t)  is  piecewise  constant  in  time.  Specific- 
R  O 

ally,  take  8  ■  0  and  let 


a(t)  -  aQ 

for  t  <  t  <  t. 
o  —  1 

a(t)  -  a  jt  a 

1  o 

for  t^  _<  t  <  “ 
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Then  Eqs.  (11)  and  (14)  and  the  Lemma  can  be  used  to  show  that 


V”'  -  “o'1  W  *  (“l'1  -  “o'1)CR(to)e 


-“oCt.-t  ) 


which  Is  Incompatible  with  the  conclusion  of  the  Theorem  as  expressed  in 
Eq.  (13)  except  when  »  ctQ. 

The  discussion  in  the  preceding  paragraph  suggests  the  conjecture 
that  in  the  general  case  represented  by  Eqs.  (9)  and  (10)  there  will  be 
an  optimal  ventilation  policy  a(t)  for  each  schedule  8(t)  of  introduction 
of  contaminant  vapor  into  R.  Determining  the  optimal  ventilation  policy 
a(t)  subject  to  appropriate  technological  constraints  on  the  feasible 
ventilation  rate  policy  may  involve  an  application  of  modern  control 
theory,  the  calculus  of  variations,  or  other  methods.  If  the  contaminant 
vapor  introduction  rate  schedule  8 (t)  is  not  known  exactly  but  must 
instead  be  estimated  using  slow  and  inaccurate  instruments,  we  conjecture 
that  there  will  still  be  an  optimal  ventilation  policy,  although  its 
determination  may  be  difficult.  The  further  development  of  these  ideas 
and  conjectures  is  beyond  the  scope  of  this  document. 
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PARTICULAR  CASES 


Suppose  that: 

1.  There  is  no  filtration,  so  E  ■  0.  (Then  to  satisfy  the  conditions 
of  the  Theorem  FQ(y)  must  be  a  positive  constant.) 

2.  The  initial  vapor  concentration  CD(t  )  ■  0. 

K  O 

3.  There  are  no  internal  contaminant  vapor  sources,  so  GQ(t)  ■  0. 

4.  The  vapor  concentration  CQ(y)  in  the  influx  stream  is  the  same 
as  in  the  outside  air. 

This  situation  corresponds  to  Example  1  in  the  Introduction.  By  Eq.  (15) 
the  total  vapor  dosage  inside  the  shelter  will  be 

dr(«)  -  f  CQ(y)  dy, 

*"0 

which  is  the  same  as  the  total  vapor  dosage  outside  the  shelter,  as 
asserted  in  the  Introduction. 

This  mathematical  result  does,  however,  neglect  some  potentially 
important  considerations.  For  example,  when  the  shelter  volume  V  is 
large  and  the  leakage  flow  is  slow,  the  vapor  concentration  inside  the 
shelter  may  be  much  lower  than  outside,  although  (since  we  know  that 
the  total  vapor  dosage  will  be  the  same  inside  as  outside)  the  inside 
vapor  concentration  must  persist  for  a  correspondingly  longer  period 
of  time.  In  reality,  however,  some  contaminant  vapor  is  adsorbed 
on  surfaces  within  the  shelter  and  the  effect  of  this  is  relatively 
more  important  the  larger  the  shelter  (because  it  will  have  a  larger 
adsorbing  surface)  and  the  lower  the  vapor  concentration  (because 
adsorbtion  will  cause  a  proportionally  larger  decrease  in  the  vapor 
concentration  levels) .  Also,  since  the  human  body  does  neutralize 
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and/or  eliminate  contaminants  at  some  rate,  sufficiently  low  vapor 
concentrations  will  not  accumulate  and  produce  intoxif ication.  In 
addition,  it  is  possible  for  the  external  vapor  concentration  to  even¬ 
tually  fall  below  the  internal  vapor  concentration.  Should  this  occur, 
shelter  occupants  should  either  abandon  it  or  ope1-,  it  up  and  air  it  out 
(an  action  that  increases  the  ventilation  rate  and  deliberately  violates 
its  assumed  constancy) . 

Next  suppose  that 

1.  The  initial  vapor  concentration  CD(t  )  =  0. 

K.  O 

2.  The  influx  vapor  concentration  Cq (y)  =  0. 

This  corresponds  to  Example  2  in  the  Introduction.  By  Eq.  (15)  the 
total  vapor  dosage  inside  the  shelter  will  be 

DR  (ot)  -  Go(y)  dy  '  [Fo  +  FAE] 
o 

i.e.,  it  is  the  ratio  of  the  total  mass  of  agent  that  eventually  vaporizes 
divided  by  the  shelter's  volume  ventilation  rate,  as  asserted  in  the 
Introduction. 

This  mathematical  result  does,  however,  neglect  some  potentially 
important  considerations.  For  example,  if  the  outside  vapor  concentra¬ 
tion  is  lower  than  that  inside,  then  shelter  personnel  should  abandon 
the  shelter  or  air  it  out  (which  deliberately  violates  the  assumed 
constancy  of  the  ventilation  rate).  Although  the  mathematical  result 
holds  for  any  size  shelter,  larger  shelters  tend  to  have  larger  volume 
ventilation  rates.  However,  they  also  tend  to  have  more  traffic  entering 
and  leaving  them  and  so  may  experience  more  liquid  or  solid  contamination. 
It  is  not  clear  which  of  these  effects  predominates.  When  the  evaporation 
rate  is  sufficiently  slow  relative  to  the  volume  ventilation  rate  the 
vapor  concentration  level  will  be  low  and  adsorption  and  biological 
resistance  can  become  important  factors. 
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CONCLUSIONS 

We  have  demonstrated  a  useful  and  interesting  general  result  in  the 
theory  of  vapor  contamination  hazards  to  collective  protection  systems. 
The  combination  of  the  generality  of  the  result  and  the  simplicity  of 
its  proof  may  be  new.  Whether  new  or  not,  the  general  result  and  its 
simple  proof  deserve  to  be  more  widely  known.  The  general  result  for  a 
constant  ventilation  rate  was  illustrated  by  showing  how  Examples  1  and 
2  are  obtained  as  particular  cases . 

We  have  also  shown  that  the  general  result  need  not  hold  when  the 
ventilation  rate  fluctuates.  Thus  there  presumably  are  optimal  ventila¬ 
tion  policies  or  schedules,  although  the  further  development  of  this 
suggestion  is  beyond  the  scope  of  this  document. 
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